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Dynamic binding of driven interfaces in coupled ultrathin ferromagnetic layers 
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We demonstrate experimentally dynamic interface binding in a system consisting of two coupled 
ferromagnetic layers. While domain walls in each layer have different velocity- field responses, for 
two broad ranges of the driving field, if, walls in the two layers are bound and move at a common 
velocity. The bound states have their own velocity-field response and arise when the isolated wall 
velocities in each layer are close, a condition which always occurs as if ^ 0. Several features of the 
bound states are reproduced using a one dimensional model, illustrating their general nature. 

PACS numbers: 75.70.Cn, 75.60.Ch, 75.78.Fg 



o 

a 
a 

O 

o 



> 

(N 

in 

o 
o 



% 



Moving elastic interfaces are a common feature of 
many non-equilibrium phenomena. These include crys- 
tal growth [1| , wetting [2|, [sl , combustion [j, |5[ , vortex 
motion in superconductors [6], mechanical fracturing [7] 
and switching in ferromagnetic and ferroelectric materi- 
als |8Ml2|. However, the majority of work examining the 
universal dynamics of these interfaces has been limited to 
studies of single interfaces, and it is only quite recently 
that the problem of coupled interfaces has begun to be 
addressed fTsHTl. 



Up until now, experimental studies have been re- 
stricted to the particular problem of two interfaces mov- 
ing through a single mediumgenerally while separated 
by a finite lateral distance [3|, ll6| . In this Letter we con- 
sider interface dynamics in a novel type of experimen- 
tal system consisting instead of two coupled, but physi- 
cally separate, media. The interfaces, magnetic domain 
walls, move through two ferromagnetically coupled ultra- 
thin ferromagnetic layers under the action of an applied 
driving field, H. The central result of this work is clear 
evidence of a dynamic binding of the domain walls in the 
two media for certain finite ranges of H. Due to the two 
layers having different disorder strengths and thicknesses 
[9] , in the absence of coupling the domain walls in each 
layer have different velocity- field responses. Experimen- 
tally however, we find that at two certain values of the 
driving field, H = and ii" = ii"* > 0, the velocities in 
each layer are the same. We show that due to the inter- 
layer coupling, a dynamic binding of the walls appears 
over finite ranges of the applied field near the crossing 
points, H = and H"" . In these field ranges, walls in 
the separate layers move together at a common velocity. 
Notably, the bound states are characterized by their own 
velocity-field response, many features of which can be re- 
produced by a one dimensional model which takes the 
strength of the interlayer coupling into account. 



Our system, shown in Fig. 1(a), consists of two ul- 
trathin weakly disordered ferromagnetic Co layers with 
perpendicular anisotropy which interact via a net fer- 
romagnetic (FM) interlayer coupling [17] of energy J. 
Sandwiched between Pt, such ultrathin Co layers are well 
established as model systems for studying the dynam- 
ics of one dimensional (ID) interfaces moving through a 
weakly disordered two dimensional (2D) medium (walls 
are narrow ~10 nm) [sMllj. The multilayer stack, hav- 
ing structure Pt(4.5 nm)/Co(0.5 nm)/Pt(3 nm)/Co(0.8 
nm)/Pt(3.5 nm) was sputtered at room temperature onto 
an in-situ etched Si/Si02 substrate. The two Co layers 
switch together during hysteresis, consistent with an FM 
coupling [18]. The magnetically 'hard' thick 0.8 nm Co 
layer has a stronger depinning field as compared to the 
'soft' 0.5 nm Co layer which, as will be seen below, re- 
sults in slower field induced domain wall propagation at 
low field. 

Domain wall velocities, v{H), with ii" applied perpen- 
dicular to the film plane, were determined quasi-statically 
from high resolution (0.4 /im) far-field polar magneto- 
optical Kerr effect (PMOKE) microscopy images. Both 
magnetic layers were first saturated in a strong negative 
field (lii'l = 1 kOe to 4 kOe). Short positive field pulses 
(~1 kOe over ~100 ns) could then be used to nucleate 
isolated 'spin up' domains in the hard layer and aligned 
'spin up' domains in both layers [Fig. 1(b)] [19]. These 
latter aligned domains allow us to study bound domain 
wall dynamics. Domain configurations were imaged both 
before and after the application of a second field pulse 
applied to drive wall motion. These two images were 
subtracted from each other [eg. Figs, l(c-h)] and the av- 
erage wall displacement and velocity were determined. 
Further details regarding this method can be found in 
Ref. [9]. 

The characteristic field dependence of the velocities of 
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FIG. 1: (a) 'Hard' (h) and 'soft' (s) Co layers, sandwiched 
between Pt layers, interacting via a ferromagnetic inter layer 
coupling, J. Domain walls separate 'spin up' and 'spin down' 
domains, (b) PMOKE image of magnetic domains in the 
hard and soft layers showing the relative orientations of the 
magnetization in each domain: large (small) arrows for the 
hard (soft) layer. 'Spin up' domains may either be isolated in 
the hard layer (mid- level gray) or exist in both layers (black). 
Subtracted PMOKE images for aligned (c-e) and hard layer 
(f-h) walls (dark regions correspond to areas swept out by the 
domain walls). Field pulse amplitudes and durations are: (c) 
112 Oe 367 s (d) 638 Oe 3x700 ns (e) 1120 Oe 250 ns (f) 155 
Oe 45 ms (g) 670 Oe 500 ns (h) +1040 Oe 2x80 ns. Black 
scale bars are 10 /xm long in images (b-h). 



isolated hard and soft domain walls, v/j, and Vg, in the 
absence of any interlayer coupling are plotted in Fig. 2(a). 
In both layers, wall motion at low field is consistent with 
a thermally activated creep regime [sl, |20| : 



v{li) = vq exp 



Uc [H, 



knT 
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1/4- 
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as demonstrated in the hiVh^s{H~^^^) plot in Fig. 2(b). 
In Eq. (1), Uc/^bT is related to the disorder-induced 
pinning energy barriers, vq is a numerical prefactor and 
1/4 is the universal dynamic exponent for a ID interface 
moving in a 2D weakly disordered medium. The slope 
of \nv [H~^^^^ is equal to —HjUc/kBT. It increases 
with the disorder strength and is higher for the thicker 
hard layer [9,] . Vh was measured directly in the hard layer 
of this film. However, since we could not nucleate isolated 
soft layer domains, Vg was measured in the soft 0.5 nm 
Co layer of a sample with a 4 nm spacer (resulting in a 
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FIG. 2: (Color online) Isolated domain wall velocities in the 
hard {vu) and soft (vs) layers in the absence of coupling plot- 
ted as (a) v{H) and (b) \nv{H-^^^). 



weak antiferromagnetic interlayer coupling [21|) but an 
otherwise identical structure to that of the film studied 
here. 

We now discuss how ferromagnetic interlayer coupling 
may theoretically give rise to bound states. Existence 
of bound states relies on two important features of our 
system. The first feature is that the Vh{H) and Vs{H) 
curves [Fig. 2(a)] cross at two field values H = H* ^ 870 
Oe and dX H = {). Note that this second crossing point is 
universal since v ^ as ilf ^ 0. The second feature, as 
illustrated in Fig. 3(a), is that the ferromagnetic inter- 
layer coupling induces an attractive interaction between 
walls in each layer. This is mediated by effective coupling 
fields, Wj = J/{tiMg) where M^ and U are the satura- 
tion magnetization and thickness of layer i = h^s [22|. 
These fields act to align the magnetization of the domains 
in each layer and hence the domain walls (eg. [23]). 

Ferromagnetic (attractive) coupling therefore com- 
petes with the tendency of the domain walls to move sep- 
arately. If Vh ^ Vs, the former prevails over the latter and 
a bound state is favored. On these very general grounds, 
from Fig. 2(a) we expect three regimes which are shown 
schematically in the left hand column of Fig. 3(b). (i) 
For H ^ and H ^ ilf*, walls are bound and propa- 
gate together at a velocity vi){H). Far from these field 
values, walls move separately, either (ii) with Vs > Vh 
or (iii) with Vh > Vs- Following Fig. 3(a), laterally sep- 
arated walls will move under the action of a total field 



[21|, |2J|, H ± Hj' , according to their relative positions: 
H -\- Hj' (increased v) for the trailing wall and H — H^ 
(decreased v) for the leading wall. 

Domain imaging was always carried out after removal 
of the driving field. Upon the removal of i7, walls which 
are separated during their field driven motion will relax 
towards each other under the action of an effective field 
equal to ±i^j as shown in the right column of Fig. 3(b). 
From the data in Fig. 2(b) (see the blue and red lines cor- 
responding respectively to \n.Vs{Hj) and \iiVh{Hj)) we 
find that Vs{Hj) ^ 1 m/s ^ lO'^^ Vh{H^), indicating that 
effectively only the soft layer wall will move significantly, 
relaxing rapidly to the hard layer wall position in a time 
much shorter than the ~10 s needed for image acqui- 




(i) 
Bound 
motion 

(ii) 
Soft layer 
wall leads 



(iii) 
Hard layer 
wall leads 



Driven dynamics (H>0) 
h| — m--^ ^ 



Post-pulse relaxation (H=0) 



hh 



hH^ 



->v,(H) 

>v,(h+h;) 



hh 



h| 

s I — o- 



► V3(H-H3) 
>v,(H-H5) 



"Vs(-Hj) 



VgCH+HJ) 



hh 



"VsC+HJ) 



FIG. 3: (a) The sign of the couphng fields, Hj' , is such that 
walls in the two layers are driven towards one another (white 
arrows). Positive applied fields will tend to drive each wall 
to the right and in the ID model the average positions of the 
walls (filled circles) are considered, (b) The total field acting 
on the walls during driven motion is if ± Hj'^ i^j'^ — for 
aligned walls). For certain field ranges neither wall can move 
faster than the other, resulting in a dynamic bound state with 
Vh — Vs — Vb (i). Outside these field ranges, unbound motion 
will occur with Vs > Vh (ii) or Vh > Vs (iii). In these latter 
cases, when the applied field is removed, the soft layer wall 
will relax to the position of the hard layer wall under the 
influence of =bif}. 



sition. This results in an apparent displacement of the 
aligned walls corresponding to that of the hard layer wall. 

In light of the above considerations we now examine 
experimental results for the coupled dynamics in the low 
field regime. In Fig. 4(a) we plot InVa as a function 
of H~^^^ where Va is the velocity of aligned walls as 
determined from their experimentally observed displace- 
ments. The crucial feature is a distinct change in slope 
dit H ^ 250 Oe. Above this field, the measured Va dy- 
namics correspond to those of hard layer walls driven 
under the action of a positive driving field and a positive 
Hj (the Vh{H + Hj) data). This is what is expected 
for unbound walls with Vs > v^^ the situation shown in 
Fig. 3(b,ii). Below this field, unique dynamics are ob- 
served with Va{H) ^ Vh{H ± Hj). Wall motion here is 
shown below to be consistent with bound dynamics. Ad- 
ditionally, since InVa ex H~'^^^^ this regime is consistent 
with a bound creep regime. Interestingly, the energy bar- 
riers for the bound walls appear to be defined by disorder 
in the hard layer (compare the slopes of the In Vh [H~^^'^) 
and \nva [H~^^^) data below 250 Oe) with the increased 
velocity for the bound walls consistent with a larger vq 
[seeEq. (1)]. 

Unbound motion above H ^ 250 Oe can be sustained 



only if Vs{H - Hj) > Vh{H + H^) [Fig. 2(b,ii)]. This 
occurs for H > Hd where 



Vs{H,,-H}) = Vh{H,,^H!}). 
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FIG. 4: (Color online) Va{H) corresponds to dynamics de- 
termined experimentally from walls aligned at remanence. 
Vh{H) and Vs{H) are isolated hard and soft layer velocities 
in the absence of coupling [Fig. 2]. Hard layer dynamics in 
the presence of a positive coupling field (vh{H + Hj)) and 
a negative coupling field {vh{H — H})) are also shown. Low 
field creep dynamics are shown in (a). Fits to the linear low 
field Vh,s,a{H) data are plotted as solid black lines. Vb(H) are 
bound wall velocities calculated using Eq. (3) with experimen- 
tally determined coupling fields. Full velocity field curves are 
shown in (b). The calculated Vb{H) for the high field bound 
state, also determined using Eq. (3) is shown as a red line. 
The vertical dotted lines correspond to the predicted lower 
and upper limits for the high field bound state. 



After fitting the low field data in Fig. 2 with Eq. (1), we 
can use Eq. (2) to estimate Hd = 254 Oe where we have 
used experimentally determined values for the coupling 



fields: H^ 



120 Oe and Hj = 220 Oe [25]. This value 



(2) 



for Hci agrees very well with the field at which a slope 
change is observed in the \nVa{H~^^^) data in Fig. 4(a) 
(see the vertical dotted line dit H = 254 Oe). For H < 254 
Oe, the coupling fields prevent either wall leading the 
other, thereby binding them [Fig. 3(b,i)]. 

The bound wall velocity can be predicted by modeling 
the coupling fields with hy^{d) = Hf^ tdinh{d/A) where 
d = ys — Vh ^^ the distance separating walls [Fig. 3(a)]. 
The out of plane coupling field is expected to be propor- 
tional to the out of plane component of the magnetiza- 
tion, m^, and this model describes m^ for a Bloch domain 



wall of width A (eg. [26|) where we use the same A for 
both layers. In essence, A simply determines the length 
over which the coupling 'force' changes sign. In place of 
Eq. (2), the condition for bound motion is then 



2007JHLPEZ). Thanks to A. Mougin for a critical read- 
ing of the manuscript. 



Vs{H - h^jid)) = Vh{H + h'}{d)) . 



(3) 



We can numerically solve Eq. (3) for d using the Vs^h{H) 
data in Fig. 2 together with H^ . If d can be determined 
for a given H^ Vh = Vs = v^ and v^ can be obtained 
simply using Eq. (3). Results are plotted in Fig. 4(a). 
Experiment agrees well with the predicted v^ for fields 
below 254 Oe where bound motion is expected. Note 
that our result does not depend on A (therefore giving 
a model with no free parameters) and that for a given 
field, no bound state exists if Eq. (3) has no solution for 
d (therefore providing an alternative method to find Hd 
which is consistent with that discussed above). 

In Fig. 4(b), a second bound state, also with a unique 
velocity-field response, is identified around H* ^ 870 Oe. 
Here, a change in slope of the Va{H) data is again ob- 
served due to a deviation away from the Vh{H ± Hj) un- 
bound behaviour. Away from this second bound regime, 
the velocity-field response of the aligned walls compares 
well to that of unbound walls with the hard layer wall ei- 
ther trailing {H < 840 Oe, Va ~ Vh{H^H^), Fig. 3(b,ii)) 
or leading {H > 1190 Oe, Va - Vh{H-Hf), Fig. 3(b,iii)). 
We again use Eq. (3) to predict this bound state's velocity 
vi,{H) and its limit fields, Hc2 ^ 600 Oe and Hc2. ~ 1050 
Oe, all shown in Fig. 4(b). While our model appears 
to capture the essential physics of the stability and dy- 
namics of the bound states, at high field it predicts a 
bound velocity which lies between the measured Vh{H) 
and Vs{H) whereas the observed bound velocity is higher. 
We note that our simple calculation of the velocities does 
not allow for any effects that may appear in a full 2D 
treatment of the problem (eg. elasticity) nor do we con- 
sider dipolar fields generated at the domain walls [27]. 

In conclusion, we report the discovery of a dynamic 
binding of driven interfaces with general characteristics 
that can be understood to be a consequence of attractive 
interactions and velocities that match at some driving 
force. Since v ^ always for a vanishing force, the 
crossing point at ilf = is universal and there may be 
analogies with other systems, for example coupled vortex 
motion in superconductors [28|, l29| and bound solitons 
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